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Abstract 

The binding energy of a negatively charged hydrogenic impurity with on- and off- 
center position in a spherical Gaussian quantum dot was calculated with the con- 
figuraction interaction method. Our calculations show that Ef, is always positive for 
on-center impurities with a maximum near to the radius for one-electron stability of 
the potential well Re- For off-center positions the binding energy can assume nega- 
tive values within a range of the quantum dot radius, thus indicating the instability 
of the system 
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1 Introduction 

Doping of semiconductor crystallites of nanometer size, or quantum dots 
(QDs), allows tuning the transport, electric, optical and magnetic proper- 
ties for the purpose of tailoring proposed quantum devices [1]. Incorporation 
of impurities into QDs provides charge carriers that strongly modifies those 
properties Neutral and negatively charged shallow donor impurities {D^ 
and D~ centers) in semiconductors are the analogue of the H atom and the 
H~ ion in atomic physics, i.e., one and two electrons bonded to a positively 
charged Coulomb center, respectively. In particular, D~ centers are the sim- 
plest system where correlation effects can play a role. The binding energy of 
a center in QDs has been studied with different confining potential shapes 
and calculation methods P][5|5f7f51l5] . Many of them assume the impurity to 
be at the center of the QD. Nevertheless, the position of the impurities 
was shown to strongly affect the binding energy p[7ll8]|9] . Other properties 
also show such a dependence; for instance, the calculated optical-absorption 
spectra of homogeneously distributed centers, show an absorption edge 
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associated with transitions involving impurities at the center of the well and a 
peak related with impurities next to the edge of the dot . Also the effect of 
parabolic confinement on the binding energy of shallow hydrogenic impurities 
in a spherical QD of a widegap semiconductor, such as GaAs, as a function 
of the impurity position for different dot sizes, was studied |7]. The binding 
energy of an off-center neutral hydrogenic impurity in a spherical quantum 
dot has been studied by using finite-depth spherical well [8] and Gaussian 
confining potentials [9]. The binding energy [10] and the energy levels of the 
ground and the excited states of spin-singlet and spin-triplet configurations 
have been calculated variationally by assuming a square finite-well confining 
potential [11] . Since the experimental demostration of the existence of built-in 
centers in doped multiple quantum wells structures [12], a number of works 
considered the binding energy of on-center negatively charged impurities, un- 
der different confining potentials [T31ll4fl5lll6lll8] . Xiew proposed a procedure 
to calculate energy spectrum of D~ centers in disk-like QDs with a parabolic 
lateral confining potential. He found that there exists a critical radius 
such that if -R < -R'^ the D~ configuration is stable p^. Pandey et al. studied 
the dependence of the binding energy of D*^ and D~ centres on the confining 
potencial shape by using the local density approximation [13]. A Gaussian 
confining potential, having finite depth and range, has been suggested as way 
to take into account effects of non-parabolicity in the QD potential for both 
one- and few-electron systems [T9ll20f21] . The energy spectra oi D~ centres in 
disk-like Gaussian quantum dots were calculated in Ref. [15]. Recently Xie cal- 
culated the binding energy of an on-center D~ donor in a Gaussian potential 
[TB] and Sahin showed that the use of the exchange and correlation potential 
is necesary, within the the local density approximation of density functional 
theory, for obtaining correct results [17] . In Ref. [18] the binding energy of an 
off-center D~ impurity in a two-dimensional parabolic QD was addressed by 
using finite-difference and fractional dimension methods. 
To our knowledge, the issue of an off-center D~ donor in a spherical QD has 
not been addressed. Therefore, the purpose of the present work is to study the 
binding energy oi a. D~ center in a spherical QD as a function of its position, 
and to explain this dependence in terms of a simple model. 



2 Theory 



We consider two electrons bonded to a shallow donor impurity in a spherical 
QD of radius R and depth Vq. The impurity is located at the position d and a 



Gaussian confining potential V{r) 



is assumed for the QD. The 



Hamiltonian, in the effective mass approximation, can be written as 
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where W{d,r) = — |r — d|~^ is the electron-donor Coulomb potential. We 
use the donor Bohr radius an = {e/m*)a-Q as the unit of length and the donor 
effective atomic unit a.u.* = (m*/e^) Hartree as the unit of energy. The binding 
energy of the D~ center is defined as [TOj 

E, = E{D'') + E{e)-E{D-), (2) 

where E[D^) is the energy of the neutral impurity in the QD, E{e) is the 
energy of an electron in the QD wihout the impurity, and E{D~) is the energy 
of the D~ in the QD. The energies E{D^) and E{e) of the one-electron systems 
are calculated by direct diagonalization. The calculation method was reported 
elsewhere [21]. The energy E{D~) of the two-electron system were calculated 
with the configuration interaction (CI) method [22], where the eigenvectors of 
two electron hamiltonian Eq. ([T]) are expanded in terms of the two-electron 
Hartree-Fock ground state and its single and doubly excitated configurations 
(Slater determinants), expanded in a single-particle Cartesian Gaussian basis 
set 

Iff =x"'y''zPexp{-air^), (3) 

where i = m + n + p is the angular momentum of the function, and the ctj 
are properly chosen exponents [21j. A basis set 4s4p4d centered in the QD, 
and a 8s7p2d basis set centered at r = d, similar to other previously used for 
describing the weakly bonded H~ ion and its polarizability, was also added 
for taking into account the donor center [23j. The total spin symmetry of 
the configurations considered were restricted to 5* = as in previous works 
[Tl|13|15l| . A potential depth of Vq = 25 a.u.* was kept throughout the present 
work. 



3 Results 

The binding energy, Eq. ([2]), calculated with the CI method is shown in Fig. 
[T] with empty circles as a function of the QD radii R. Four impurity positions 
were considered, namely, d = 0,0.3,0.6 and 1.0 ao. The results show that 
for on-center position [d = 0) Eb is always postive with a maximum nearly 
i?c — 0.2a/3- This maximum binding energy at this critical radius Rc, is related 
to the fact that for every Vq, there is a minimum R^ where an electron can 
be stable in the QD [2T] . This result is in qualitative agreement with previous 
works that treated on-center D~ donors [T0|14fl7j . It should be mentioned, 
however, that they differ from a recent calculation by Xie [T6]. At small dis- 
tances from the potential center {d = 0.3a£)), the maximum is less pronounced 
and there is a minimum at Rc- For larger values of d (0.6 ajy and lao), there 
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still exists a minimum at Rc, such that the larger d, the more negative the 
minimum becomes while the maximum becomes flatter. Also the binding be- 
comes negative for radii R ^ d, and positive for R> d. Hencefore, the larger 
d, the wider the radii range where the binding energy is negative. We also 
preformed Hartree-Fock calculations, not reported here, whose results show a 
similar trend. The correlation energies were found in the range of —0.032 a.u.* 
for R = 0.25aD, to —0.037 a.u.* for R = lOao, and weakly dependent on the 
impurity position. 

The results can be rationalized as follows. For a flxed potential depth and very 
small radius the effect of the potential becomes negligible because it cannot 
bind electrons. So, the two electrons are kept bonded due to the impurity 
Coulomb potential forming a H~ ion. The same happens for very large radius, 
where the bottom of Gaussian potential becomes flat and contributtes approx- 
imately with a constant potential — Vq. Then, Ei,{D^) Ei,(B.^) = 0.0277 a.u. 
for both i? — ^ and R oo. 

For intermediate radius {Rc ^ R ^ d), where the dot can allocate electrons, 
the impurity is outside the dot and the system could become instable. For 
very large radius {R ^ d), the system behaves like an on-center impurity, 
thus having a positive binding energy. 

A more quantitative explanation of the results can be obtained by using a 
variational estimate as follows. Consider a normalized s-type Gaussian trial 
function v^s(r) = (2a;/7r)'^/^ exp(— ar^), centered in the QD center. The en- 
ergy of the two-electron D~ center can be obtained as the expectation value 
of the Hartree-Fock Hamiltonian in the spin-singlet trial state ip{ri,r2) = 
ipsiri)ips{r2), thus giving 



where the terms within brackets are the expectation value of the kinetic energy 
Tq,, conflning potential Va and the impurity potential Wa for each electron. The 
last term is the Coulomb repulsion J between the Gaussian charge densities 
of each electron. The optimal exponent a is obtained by minimization of Eq. 
(jl]). In this way, using the optimal a, the ground state energy E{D^) can be 
estimated as E{D-) = 2{Ta+Va+Wa) + J. Analogously, E{D°) ^ T^+Va+Wa 
and E{e~) ^ Tq, + V^. Hence, the binding energy is approximately given by 
Eb{D-) = -Wa - J, that is. 



Eq. ([5]) implies that Eb{D ) depends directly on the interplay between the 
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nuclear attraction and electron-electron interaction. For d (the limit of on- 
center impurity), erf (a;) ^ 2x/y/7i and Eq. ([5]) becomes Eb = 2{y/2 — l)^a/n, 
thus showing that the on-center binding energy is always positive, in agreement 
with the CI results presented here. On the other hand, for a fixed R, as d 
increases, eTf{\/2ad)/d — > 0, and J becomes dominant, thus giving Eh < 0. 
The systems for which Eb < are not stable and is similar to a molecular 
dissociation process ending up with one electron in a QD and the other in the 
D~ . The results calculated with Eq. ([5]) are shown for comparison in Fig. [1] 
with continuous lines. As can be seen, all qualitative features of the CI curves 
are well reproduced with this simple model. It is interesting to point out that 
the use of the Hartree-Fock Hamiltonian gives an electron-electron interaction 
2J — K, where K is the exchange energy such that K = J for the doubly 
occupied ground state. Thus, Eq. (jlj) takes into account the exchange energy 
correctly. Disregarding the exchange energy would imply to add a factor of two 
to the last term of Eq. (jlj), and the binding energy for an on-center impurity 
would give Eh{D^) ~ —Wa — 2J ^ 2{\f2 — 2)y«/^ < 0, in agreement with 
Ref . [TT] . Eq. (jSj) was also used in Fig. |2] to show the change in the binding 
energy as the impurity moves from the center of the QD up to = la/). 



4 Concluding remarks 

In summary, we have calculated the binding energy of a negatively charged 
impurity with on- and off-center position in a spherical gaussian quantum dot 
with the configuraction interaction method. Our calculations show that E^, is 
always positive for on-center impurities with a maximum near to the radius 
for one-electron stability of the potential well Rc- As the impurity is displaced 
off center, the maximum of E^ decreases and a minimum near to Rc appears. 
For sufficiently large d^ Ef, assumes negative values indicating the instability 
of the system. Our results could be useful for understanding how the binding 
energy is affected by the breaking of the spherical symmetry of the potential 
well due to doping in low-dimensional systems. 
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Fig. 1. Binding energy of the Gaussian quantum dot with a negatively charged donor 
impurity at a distance d from the center for d = 0,0.3,0.6 and lao- Energies are 
given in effective atomic units and distances in donor effective Bohr radius an- The 
empty circles represent configuration interaction calculations. The continuous lines 
are results of the variational model, Eq. ([5]) . 

d=0.0 an d=0.3 a, 
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Fig. 2. Binding energy of the impurity as a function of the radius of the quantum 
dot calculated with Eq. dSJ for the impurity positions d = 0, 0.1, 0.2, 0.3, 0.4, 0.6 
and 1 a_D- Energies are given in effective atomic units. 
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